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ABSTRACT 

We present 0.9—2.5/im spectroscopy with R^800 and 1.12-1.22/tm spec¬ 
troscopy with R~5800 for the M dwarfs G1 229A and LHS 102A, and for the 
L dwarf LHS 102B. We also report IZJHKL' photometry for both compo¬ 
nents of the LHS 102 system, and L' photometry for G1 229A. The data are 
combined with previously published spectroscopy and photometry to produce 
flux distributions for each component of the kinematically old disk M/L-dwarf 
binary system LHS 102 and the kinematically young disk M/T-dwarf binary 
system Gliese 229. The data are analyzed using synthetic spectra generated by 
the latest “AMES-dusty” and “AMES-cond” models by Allard & Hauschildt. 
Although the models are not able to reproduce the overall slope of the infrared 
flux distribution of the L dwarf, most likely due to the treatment of dust in 
the photosphere, the data for the M dwarfs and the T dwarf are well matched. 
We find that the G1 229 system is metal-poor despite having kinematics of the 
young disk, and that the LHS 102 system has solar metallicity. The observed 
luminosities and derived temperatures and gravities are consistent with evo¬ 
lutionary model predictions if the G1 229 system is very young (age ^ 30 Myr) 
with masses (A,B) of (0.38,^ 0.007)Mq, and the LHS 102 system is older, 
aged 1—10 Gyr with masses (A,B) of (0.19,0.07)Mq. 
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1 INTRODUCTION 


In the last few years there have been dramatic developments in the study of very low mass 
stars and brown dwarfs — where the latter are dehned to be objects with mass below that 
required for stable hydrogen burning. Sky surveys have found a signihcant population of low 
mass dwarfs cooler than M-dwarfs. The hrst population to be identihed was the L-dwarf 
class, distinguished from the M-dwarfs by weakening VO and TiO absorption features in 
the red, and by stronger alkaline and H 2 O absorption features in the red and near-infrared 
(e.g. IDelfosse et al. IQFtI ; [Kirkpatrick et al. 1999| ; [Martin et al. 1999|) . The L-dwarfs cover 
an effective temperature range of about 2200 K to 1400 K (e.g. [Leggett et al. 200^) and 
the prototype of this class is GD 165B, discovered by [Becklin fc Zuckerman (1988|) as a 
red companion to a hot white dwarf. Several objects even cooler than the L-dwarfs were 
discovered in 1999 by the Sloan Digital Sky Survey and the 2 Micron All-Sky Survey ([Strauss 
et al. 199^ [Burgasser et al. 1999|) . They are very similar to Gliese 229B, the extremely low- 
mass companion to an M-dwarf found by [Nakajima et al. (1995[) , and are distinguished by 
the presence of GH4 absorption in the near-infrared H and K bands. These objects, known 
as T-dwarfs, have Tes ~1300—800 K and tens of T-dwarfs are now known (see [Geballe 
et al. 200^ ; [Burgasser et al. 2002[ , and references therein). 

In this paper we present an analysis of the low-mass binary systems LHS 102 and Gliese 
(Gl) 229. The primaries of these two systems are M-dwarfs which are in the proper mo¬ 
tion catalogues by [Luyten (1979[) and [Gliese fc Jahreiss (1991[ ). The secondaries were only 
recently detected. Gl 229B, the prototypical T-dwarf, was discussed above; LHS 102B is 
an L dwarf discovered by [Goldman et al. (1999[) . We present new UKIRT infrared photom¬ 
etry and spectroscopy of these objects in §2, which are analyzed using the latest models 
by Allard & Hauschildt (pdauschildt, Allard fc Baron 1999[) , described in §3. The results of 
the model comparison are given in §4, the implications for age and mass of each compo¬ 
nent are discussed in §5 and our conclusions given in §6. These two systems span the newly 
dehned low-mass spectral classes and, being binaries at known distances where each com¬ 
ponent has presumably the same chemical composition and age, are potentially very useful 
for constraining atmospheric and evolutionary models of ultracool dwarfs. 


2 OBSERVATIONAL RESULTS 
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2.1 The Sample 


Table 1 gives the LHS and Gliese or Gliese-Jahreiss catalogue numbers (|Luyten 1979| ; [Gliese 
fc Jahreiss 199 Ij) for the two systems discussed in this work. Although both primaries are 
in both catalogues we follow tradition and refer to one as LHS 102 and the other as G1 229. 
The coordinates in the Table are based on our observations at UKIRT and are accurate to 
about an arcsecond. 

Spectral types also are listed in Table 1. For the M dwarfs LHS 102A and G1 229A these 
have been taken from jMartm et al. (1999|) and [Leggett (1992|) , respectively. The T subtype 
for G1 229B is taken from the provisional classihcation scheme of [Geballe et ah (2002[) . For 
LHS 102B the indices presented by Geballe et al. imply L4.0 based on the red spectrum, 
L4.0 based on the K spectrum, and a later type of L6.5 based on the H-band water index. 
As the last region is slightly noisier than the others we adopt a type of L4.5±0.5. This is 
consistent with previously published classihcations based on red spectra only: [Martin et al. 
( 1999[ ) classify this object as L4, [Kirkpatrick et al. (2000|) classify it L5. 


Table 1 also gives the kinematic population for each system, based on the space motion 
of the primaries, where distances were obtained from the Hipparcos Gatalogue for G1 229 
( [Ferryman et al. 1999[) and from the Yale Gatalogue for LHS 102 ( [van Altena et al. 19941) . 
The population for the G1 229 system is from [Leggett (1992[) ; we determine the population 
for the LHS 102 system to be older, based on the specihcations given by Leggett for UVW 
space motions. The space motions for LHS 102A were calculated using a radial velocity from 
Rodgers fc Eggen (1974[) ; a more recent but quite different value for the secondary from 
Basri et al. (2000[) gives the same end result — that the V motion at —73 (for LHS 102A) 
or —77km/s (for LHS 102B) is too large for the “young” disk. 


2.2 New Photometry and Observed Colours 

New photometry has been obtained for three of the stars in the sample. In 1999 August, 
IZJHK was obtained for the LHS 102 system using the UK Infrared Telescope (UKIRT) 
with the UFTI camera, a 1024x1024 imager with O'.'OO pixels. The I and Z Liters are non¬ 


standard and have half-power bandwidths of 0.72—0.93 pm, and 0.85—1.055 pm. [Landolt 
( 1992[) standards were used to calibrate the / data and provisional Sloan standards were 


used for the Z data ( [Krisciunas et ah 1998[) . For both I and Z there are strong colour terms 
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between the UFTI system and the Landolt and Sloan systems; the results are given in Table 
2 in the natural UFTI system. 

The UFTI JHK hlters are part of Manna Kea consortium hlter set ([Simons fc Tokunaga 
20021 ; [Tokunaga fc Simons 2002]) . UKIRT faint standards ([Hawarden et ah 2001[) were used 
to calibrate the photometry. There are colour differences between the IRCAM3 system used 
to establish the UKIRT faint standards and the natural UFTI system, dehned by the new 
hlter set (for more discussion see [Hawarden et al. 2001[ ; [Leggett et al. 200^) . In Table 2 we 
list JHK in the natural UFTI or Mauna Kea consortium (also known as MKO) system. 

New L' photometry was obtained for LHS 102 A and B, and G1 229A, in 1999 September, 
using UKIRT’s IRCAM/TUFTI camera, a 256x256 imager with O'.'OS pixels. Bright UKIRT 
standards were used to calibrate the photometry. The L' magnitudes of late L dwarfs and T 
dwarfs are different for the Mauna Kea consortium hlter and the hlter previously installed in 
IRC AM and referred to as the UKIRT hlter (see [Leggett et al. 200'2| [Stephens et al. 20011 ); 
the results in Table 2 are on the MKO system. 

Table 3 gives for both systems the distances taken from the parallax catalogues ([Per¬ 
ryman et al. 199^ ; [van Altena et al. 1994[ ) in the form of distance moduli, and tabulates 
IZJHKL' colours for the sample. In this Table we give I magnitudes on the Cousins system. 
The value for Cl 229A is taken from the compilation by [Leggett (1992[) . For the other three 
objects the Ic value is synthesized from the hux-calibrated spectral energy distributions 
using a hlter prohle from [Bessell (1990[) (red spectra from the literature have been added 
to the infrared spectra, see §2.3). For Cl 229B we extended the available spectroscopic data 
shortwards of 0.8/rm to 0.7pm using as a template the energy distributions for the similar 
object in [Strauss et al. (1999[) . The uncertainties in I for Cl 229B and LHS 102A,B are 
~ 10%. We suggest that these values for I are preferable to those reported by [Coldmaii 
et al. (1999[) for LHS 102 A and B, which are inconsistent with our hux distributions, as 
is the value for LHS 102A reported by [Rodgers fc Eggen (1974[) ; i.e. adoption of either of 
these earlier results would lead to a discontinuity in the observed red to near-infrared hux 
distribution. The Z magnitudes for Cl 229 A and B are synthesized from their spectral 
energy distributions (the primary would saturate UFTI and the secondary is contaminated 
by the primary in imaging mode); the uncertainties are ~ 10%. The MKO system JHK for 
Cl 229A are synthesized from our hux-calibrated spectral energy distributions (the spectra 
were hux calibrated using the photometry reported in [Leggett 1992[) ; the uncertainties are 
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2.3 New Spectroscopy and Spectroscopic Sequences 


New R~5800 spectroscopic data for 1.12—1.22/nn for G1 229A and LHS 102A,B were ob¬ 
tained during 1999 August—November, using the CGS4 spectrometer on UKIRT. Telluric 
features were removed by ratioing with a nearby A- or F-type star (after removing hydrogen 
lines present in the ratio star spectrum). The data were obtained in two spectral segments 
which were then joined around 1.16-1.17 /im. Data with R~3000 over narrow intervals in 
the J, H and K bands for G1 229B were taken from Saumon et ah (2000) and also used in 
this analysis. 

Spectra at R~800 covering 0.86—2.52 /rm for G1 229A and LHS 102A, and 0.83—2.52 jam. 
for LHS 102B, were obtained on 1999 September and August, also using GGS4 on UKIRT. 
Telluric features were removed by ratioing with nearby A- or F-type stars (after removing 
hydrogen lines present in the ratio star spectrum). The spectra consisted of several segments, 
each segment was flux calibrated using our photometry and conjoined with adjacent segments 
to form spectral energy distributions. We excluded the range 1.36—1.41 jim where terrestrial 
H 2 O absorption makes the data too noisy to be useful. 

The 0.84—4.15 jim spectra for G1 229B is taken from [Leggett et ah (1999|) who re¬ 
flux-calibrated the red and L-band spectra of Pppenheimer et ah (1998|) and the near- 
infrared spectra of [Geballe et ah (1996|) using photometry. These data have been further 
rehned by putting the Oppenheimer et ah data onto vacuum wavelengths and calibrating 
them using the HST photometry of [Golimowski et ah (199^ ). We also use the R~480 4.5— 
5.1 jim spectrum of G1 229B from [Noll et al. (1997]) . This spectrum is not flux calibrated 
photometrically, and the flux level may be uncertain by ~30% (due to possible differences 
in slit losses between G1 229B and the calibration star). 

Spectra shortwards of 0.86 jim have been taken from [Henry et al. (199^ ) for G1 229A 
(R^idOO), from [Henry et al. (2002[) for LHS 102A (R?^700), and from [Martm et al. ( 1999[) 
for LHS 102B (Rf«700). The blue ends of the GGS4 spectra were calibrated using our Z 
and J photometry and the red spectra were scaled to match the GGS4 data in the region of 
overlap. 


2.4 Integrated Fluxes and Bolometric Corrections 

Table 3 gives integrated fluxes for the sample, expressed as fluxes at the Earth, bolometric 
magnitudes and intrinsic stellar luminosities. The updated results for G1 229B are expressed 
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on the MKO system (also presented in [Leggett et al. 20021) . For the remaining objects, the 
integrated fluxes were obtained by integrating the spectroscopic data over wavelength, and 
adding the flux contributions at shorter and longer wavelengths. The shorter wavelength 
flux contribution was adopted to be a simple linear extrapolation to zero flux at zero wave¬ 
length for LHS 102B, but for G1 229A the B magnitude ( [Leggett 1992|) , and for LHS 102A 
the B and V magnitudes ( [Rodgers fc Eggen 1974[) , were used to estimate the shorter wave¬ 
length flux distributions, using an effective wavelength approach. The flux contributions at 
wavelengths beyond 2.4/rm were calculated by deriving the fluxes at L' using an effective 
wavelength approach, including the contributions between the long wavelength end of the 
K-band spectrum to the L' wavelength, determined by linear interpolation, and assuming 
Rayleigh-Jeans tails beyond LL The last is a reasonable assumption for M- and mid-L- 
dwarfs, but not for methane dwarfs, as discussed in [Leggett et al. (19^ . The uncertainties 
are 5%, 0.05 mag and 0.02dex in total flux, bolometric correction and logioL/L q, respec¬ 
tively. The luminosities are used later in our analysis, together with the derived atmospheric 
parameters, to constrain masses and ages for the binaries. 


3 MODELS AND COMPARISON TO SYNTHETIC SPECTRA 
3.1 Models 


The models used for this work were calculated as described in [Allard et al. (20dl[) , and are 
based on the Ames H 2 O and TiO line lists by [Partridge fc Schwenke (1997[) and Schwenke 


(1998[) . We stress that large uncertainties persist in these opacities for the temperature range 


of this work (see [Allard et al. 2000[) . Test calculations performed by replacing the Ames water 
with the water line list from the SCAN database ([Jprgensen et al. 2001j) , for somewhat lower 
temperatures, show relatively small effects at low spectral resolution; therefore we retained 
the Ames H 2 O line list for this analysis. The models and their comparisons to earlier versions 
are the subject of a separate publication ( [Allard et al. 200T[ ) and we thus do not repeat their 
detailed description here. 

The models have been upgraded with (i) the replacement of the JOLA (Just Overlapping 
Line Approximation) opacities for FeH, VO and CrH by opacity sampling line data from 
Phillips fc Davis (1993|) and R. Freedman (NASA-Ames, private communication), and (ii) 
the extension of our database of dust grain opacities from 10 to 40 species. In the following, 
these models will be referred to as “AMES-dusty” for models in which the dust particles 
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stay in the layers in which they have formed, and “AMES-cond” for models in which the 
dust particles have sunk below the atmospheric layer where they originally formed (and 
below the photosphere). The term “AMES grid” refers to either of the “AMES-dusty” or 
“AMES-cond” models. 

In order to investigate the effects of mixing length on the model structure and synthetic 
spectra, the AMES grid was extended to include mixing lengths £ of 1.5 and 2.0 times the 
local pressure scale height Hp for the parameter range of interest. These additional models 
were otherwise constructed with the same setup as those discussed in [Allard et ah (200ID to 
avoid systematic errors in the analysis. 

The primaries (A components) analyzed in this paper have effective temperatures that 
are too high for dust to form and to affect their atmospheres. These relatively high effective 
temperatures allow us to compare the data with the NextGen grid of model atmospheres 
for dwarfs (jHauschildt, Allard fc: Baron 1999|) and giants (jHauschildt et ah 1999|) and thus 
to compare the results obtained using older and newer model grids. The T dwarf G1 229B 
is in the regime where dust has formed and settled below the atmosphere (cf. [Adard et al 


2001|) whereas the L dwarf LHS 102B is in the regime where dust formation and opacities 
have to be considered. For each of these cooler objects we use the newer and appropriate 
models for the analysis (see §4.2 and §4.4). 


3.2 Comparison Process 

Several model atmosphere grids were generated, covering the range 500 K < Teg < 4000 K, 
3.5 < log( 5 f) < 5.5 and [M/H]= —1.0, —0.7, —0.5 and 0.0 for mixing lengths of 1.0, 1.5 and 
2.0, for a total of ~ 300 model atmospheres. Synthetic spectra generated by these models 
were compared to the observed spectra using an IDL program. In addition, the observed 
spectra were compared against the NextGen grid ( jHauschildt, Allard fc Baron 1999| ) over a 
wider range in metallicity but for only one value of the mixing length, in order to assess the 
stability of the derived parameters and to search for systematic problems. 

First, the resolution of the synthetic spectra was degraded to that of each observed 
spectrum by convolution with a Gaussian of the appropriate width, and the spectra were 
normalized to unit area for scaling. Next, for each observed spectrum the program calculated 
a quality function q, similar to a value, for the comparison with all synthetic spectra in the 
grid. The quality function is calculated by first scaling the model spectrum to the observed 
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fluxes, then mapping the synthetic spectrum (reduced to the resolution of the observed data) 
onto the grid of observed wavelength points, and, hnally, calculating 



with Wi = 0.5{f°^l + where is the (mapped) flux of the model 

spectrum, is the observed flux, and A°'^® the observed wavelength. 


We selected the models that resulted in the 3—10 lowest q values as the most probable 


parameter range for each individual star. The “best” value was chosen by visual inspection. 
This procedure allows a rough estimate of the uncertainty in the stellar parameters. Note 
that it does not eliminate systematic errors in the stellar parameters due to missing, incorrect 
or incomplete opacity sources. Table 4 lists model parameters and the goodness of £t value 
for each object, for models with a range of parameters selected to cover the probable range. 
Only the relative value of q is important, the actual value is arbitrary as it depends on the 
number of data points used and on the absolute flux levels; a decreasing value indicates a 
better £t. 

4 RESULTS 
4.1 LHS 102A 

The spectrum of LHS 102A together with the best fitting synthetic spectra generated from 
the Ames grid models is presented in Fig. |l]. The best htting models have T^s ~ 3200 K and 
log( 5 ') ~ 5.0, with solar abundances. The formal uncertainty is about 100 K in Tgs and 0.5 
in log( 5 '). Metal poor models gave a signihcantly worse fit as indicated in Table 4. NextGen 
model hts to LHS 102A result in derived parameters of Teg ~ 3400 K, log( 5 ') ~ 5.5 and solar 
abundances. However although the NextGen models include both metal-rich and metal-poor 
models, the mixing length is hxed at a = 1.0. Gonsidering the results from all sets of hts, 
we conclude that Tes ~ 3200 K, log( 5 f) ~ 5.0, solar abundances, and a = 2.0 are the likely 
parameters for LHS 102A. 

The hts to the high-resolution spectrum of LHS 102A are shown in Figs. and We have 
used the Ames grid models to generate high-resolution synthetic spectra in the parameter 
range (Teg, log( 5 f)) found from the low-resolution htting and selected the three models that 
best ht the high-resolution spectra. Many of the observed lines are reproduced quite well 
in the synthetic spectra, indicating that the derived parameters and their ranges are good. 
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There are no strong indications for non-solar patterns that can be established with conhdence 
from the data, but spectral lines from more elements are required to conhrm this result. 


4.2 LHS 102B 


In Fig. ^ we show the observed spectrum and best hts to LHS 102B using “AMES-dusty” 
models from |Ahard et ah (2001|) ; for comparison we show the results obtained using “AMES- 
cond” models (in which the dust grains are assumed to have settled completely below the 
atmosphere) in Fig. As to be expected, the “AMES-dusty” models result in substantially 
better hts to the water bands than the “AMES-cond” models, although neither reproduces 
the heights of both the J-band and K-band hux peaks. LHS 102B’s derived atmospheric 
parameters are Teg = 1900 K ± 100 K, log( 5 f) ~ 6.0, solar abundances; at this temperature 
dust formation and opacities are important ([Allard et al. 20011) . The parameters derived by 
the htting procedure using “AMES-dusty” and “AMES-cond” models are quite similar. This 
increases our conhdence in the derived set of parameters, as the “AMES-dusty” and “AMES- 
cond” models are extreme cases and partial settling models (under construction) will fall in 
between these limits (see [Allard et al. 2001] , for more details). The model parameters listed 
in Table 4 are for the “AMES-dusty” model grid. 

The “AMES-dusty” ht to the high resolution spectrum of LHS 102B is shown in Fig. |^. 
The synthetic spectra ht the atomic line prohles quite well and many molecular bands 
(hundreds of overlapping line for multiple systems) are reproduced reasonably well. As in 
the case of LHS 102A, the spectrum is well ht with solar abundances. 


4.3 G1 229A 

Fig.0 shows the best ht to the low resolution spectrum of G1 229A. The automatic htting 
procedure yielded a best htting model from the Ames grid with the parameters Teg = 3700 K, 
log( 5 ') = 4.0, [M/H] = —0.5 and a mixing length parameter a = £/Hp = 2.0. However, a 
model with a = 1.5 hts only marginally worse (Fig. for the same ehective temperature 
and gravity. The third best ht is signihcantly worse in quality than the hrst two and gives 
parameters of Tes = 3600 K, log( 5 f) = 4.5, [M/H] = —0.5 and a = 1.0. Table 4 shows that 
solar metallicity models gave signihcantly worse hts. Repeating the htting procedure with 
NextGen models gave best hts in the range of 3500 K < Tgfj < 3600 K, 3.5 < log( 5 f) < 4.0, 
—1.5 < [M/H] < —1.0. The mixing length of the NextGen grid is hxed at a = 1.0 and the 
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quality of the NextGen fits is noticeably worse, the AMES grid models (even at a = 1.0) 
£t the observed spectra much better. We conclude that G1 229A is somewhat metal poor, 
[M/H] —0.5, and that the best htting mixing length parameter a is between 1.5 and 

2.0. The latter is consistent with 3D hydrodynamical models of objects in this effective 
temperature regime ( [Ludwig et al. 200^) . 

The comparison of the Ames grid models to the high resolution spectra is shown in 
Figs. 1^ and ^ Fig. indicates that scaled solar metallicities may not be a good assumption 
for G1 229A: the lines of Fe I are systematically too strong whereas the lines of Na I are too 
weak in the synthetic spectra. The K I A1177 doublet is too weak in the models (Fig. but 
this comparison is uncertain because just longward of the K I doublet there is some blending 
with Fe I lines. We note that jSchiavon et al. ( 1997] ) derived [Fe/H]= —0.2 for G1 229A based 
on the FeH band at l/im, however [Mould (1978|) found [M/H] = +0.15 based on Fourier 
analysis of Al, Ga and Mg infrared lines. These conflicting results in the literature may 
support a non-solar metallicity pattern for this star. 


4.4 G1 229B 

The best hts for G1 229B are shown in Fig. For this object all model htting was done 
by hand as the strong water and methane absorption bands gave rise to erroneous results 
using the automated routine. The range of parameters for G1 229B is Teg = 1000 K ± 100 K, 
log(fi') 3-5 and [M/H] —0.5, putting it hrmly into the “AMES-cond” regime; “AMES- 

dusty” models cannot fit its spectrum at all. The results corroborate the hts of G1 229A with 
lower metallicity models so that the G1 229 system does indeed seem to have low metallicities 
despite it being a member of the kinematic young disk. Solar abundance models do in general 
deliver worse hts than those with [M/H] = —0.5 (see Fig. |^, and models with even lower 
metallicities did not give better hts. Other authors have also suggested that G1 229B is 
metal-poor; [Grihith fc Yelle (2000[ ) hnd [0/H]= —0.5 based on their model analysis of the 
0.8—1.0 /im spectrum, and ^aumon et al. (200^) , in an analysis of the same infrared spectra 
studied here, hnd atmospheric parameters Tes = 870 — 1030 K, log( 5 f) = 4.5 — 5.5 and 

[M/H] = —0.5-0.1. Although young metal-poor stars are uncommon they do exist, as 

shown by [Feltzing et al. (2001 ), in particular their Figures 13(e) and 18(a). 

The hts to high resolution spectra of G1 229B are shown in Figs. pT|-[T^. For the most 
part the synthetic spectra ht the atomic line profiles quite well, but signihcant discrepancies 
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remain, in particular in the 1.5-1.7 /xm range (Fig. 0 )- These problems are likely caused in 
part by incomplete molecular line data, in particular water vapor and methane. The analysis 
roughly conhrms the results from the low resolution spectra. 


5 COMPARISON TO EVOLUTIONARY MODELS 


Table 5 lists the best £t model parameters, together with the observed luminosities of each 
of our targets. These results can be used to constrain mass and age for each component 
through evolutionary models. We used recent evolutionary models but were obliged to use 
three different models for the comparison as each covered a limited range in temperature and 
mass. For the M dwarfs LHS 102A and G1 229A we used the solar-metallicity and metal-poor 
models of [Baraffe et ah (1998| ), for the L dwarf LHS 102B those by [Chabrier et ah (2000|) , 
and for the T dwarf G1 229B those by [Burrows et ah (1997| , pOOlj) . The hrst two models 
use earlier generation “AMES-dusty” atmospheres, the last a “cond”-type atmosphere with 
the dust not contributing to the atmospheric opacity. The T dwarf comparison is based on 
a solar metallicity evolutionary model and is therefore more uncertain than the others. We 
constrain the age of G1 229B to be that of G1 229A, and use the evolutionary models to 
determine mass based on age, luminosity, temperature and gravity. These give self-consistent 
results, however Purrows et ah ( 2001 |) show that for a given luminosity or temperature, and 
hxed age, a more metal-poor brown dwarf has a higher mass than a solar metallicity brown 
dwarf, and so we regard the value derived here for the mass of G1 229B to be a lower limit. 

To determine the possible ranges in mass and age for each component of the two binaries, 
we adopted uncertainties in effective temperature of 100 K, in log( 5 ') of 0.5 dex, and in 
logioL/L 0 of 0.02 dex. The metallicity was considered to be hxed as evolutionary models 
are only available with either [M/H] =0 or [M/H] = —0.5. Assuming that the components 
of each binary have the same age, we hnd that the G1 229 system is very young (aged 16— 
45 Myr) with masses for (A,B) of (0.30— 0.45,;^0.007)Mq, and the LHS 102 system is older, 
aged 1—10 Gyr with masses (0.19,0.07)Mq. The ages are consistent with the kinematically 
implied disk populations (Table 1), although [Nakajima et al. (19951) suggests that the lack 
of Ha emission in G1 229A implies an age older than 100 Myr. We note however that the 
star is a hare star, which supports a young age. The young age of the G1 229 system gives 
rise to a lower mass for the T dwarf than previously assigned, only about 7 Jupiter masses, 
although as mentioned above this is a lower limit. The B component of the LHS 102 system 
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is most likely a brown dwarf, i.e. just on the substellar side of the the stellar/substellar mass 
boundary, based on its luminosity (see for example Figure 1 of Purrows et ah 2001|) . 


6 CONCLUSIONS 


In this paper we have presented the results of the analysis of two binary systems, one (LHS 
102) consisting of an M/L dwarf pair and one (G1 229) being an M/T dwarf system. The M 
dwarf primaries show that the best mixing length for effective temperatures above 3000 K 
is somewhere around 1.5 to 2.0, which is consistent with recent hydrodynamical models 
( Pudwig et al. 200^ . 

The LHS 102 system is best £t with solar abundance models. LHS 102B is an L dwarf 
with an effective temperature around 1900 K that is best £t with “AMES-dusty” models 
that include the effects of dust formation and opacity on the atmospheric structure and the 
emitted spectrum. More work is required on the detailed treatment of photospheric dust to 
produce a good match to the entire observed flux distribution; such work is in progress. LHS 
102B is found to be just below the stellar mass limit. 

The results for G1 229 indicate that this system is metal dehcient with [M/H] ps —0.5. 
We have determined an age for the G1 229 system of about 30 Myr, constrained primarily 
by the observed luminosity and derived effective temperature of the A component. This 
young age is consistent with the low surface gravity derived for G1 229B (see for example 
Figure 9 of purrows et al. 1997|) , and translates to a very low mass for this T dwarf of 7 
Jupiter masses. A better mass determination requires metal-poor evolutionary models of 
brown dwarfs, which are not currently available. 
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Table 1. Sample Description 


Name 

Other 

Names 

RA/Dec 

equinox 2000/epoch 1999 

Spectral 

Type 

Kinematic 

Population 

LHS 102A 

GJ 1001 

00:04:36.37 -40:44:02.5 

dM3.5 

old disk 

LHS 102B 


00:04:34.83 -40:44:06.0 

dL4.5 


Gliese 229A 

LHS 1827 

06:10:34.69 -21:51:48.9 

dMl 

young disk 

Gliese 229B 


06:10:34.86 -21:51:56.3 

dT6 


Table 2. New Photometry 





Name 

I (error) Z (error) 

(UFTI) 

J (error) 

H (error) K (error) 

(MKO) 

L' (error) 

LHS 102A 
LHS 102B 

9.55 (0.03) 
15.89 (0.04) 

9.12 (0.03) 
14.67 (0.03) 

8.62 (0.03) 
13.06 (0.03) 

8.02 (0.03) 
12.14 (0.03) 

7.73 (0.03) 
11.36 (0.03) 

7.53 (0.05) 
10.41 (0.05) 


G1 229A ... ... ... ... ... 4.06 (0.05) 


Table 3. Colours and Fluxes for the Sample 


Name 

M — m 

Ic 

Ic-Z 

Ic-J 

J-H 

H - K 
(MKO) 

K 

K-L' 

Flux 

W/m2 

“Li 

BC^ 

log 10 
L/Lq 

LHS 102A 

0.0997 

9.86 

0.74 

1.24 

0.60 

0.29 

7.73 

0.20 

1.69e-12 

10.45 

2.72 

-2.32 

LHS 102B 

0.0997 

16.68 

2.01 

3.63 

0.92 

0.78 

11.36 

0.95 

3.12e-14 

14.79 

3.43 

-4.05 

G1 229A 

1.1926 

6.11 

0.43 

1.05 

0.65 

0.25 

4.15 

0.09 

4.98e-ll 

6.78 

2.63 

-1.29 

G1 229B 

1.1926 

20.02 

3.84 

6.01 

-0.35 

0.00 

14.36 

2.14 

5.97e-15 

16.58 

2.22 

-5.21 


Notes 

Adopting Lq = 3.86e26 W and Mbol,© = 4.75 then: rrihol = —2.5 x log\Q{flux) — 18.978 
*’■. BCk = K 

Adopting Lq = 3.86e26 W, and if tt is parallax in arcseconds then: logioLl Lq = logio{flux) — 2 x logioir + 7.491 


Table 4. Atmospheric Model Fit Parameters 


Name 

reff/log(g)/[m/H] 

Quality 

LHS 102A 

3200/5.0/0.0 

7.92 

LHS 102A 

3200/4.5/0.0 

7.94 

LHS 102A 

3200/5.5/0.0 

7.98 

LHS 102A 

3200/5.0/-0.5 

8.12 

LHS 102A 

3000/5.0/-0.7 

8.29 

LHS 102B 

1900/6.0/0.0 

15.68 

LHS 102B 

1900/5.5/0.0 

16.21 

LHS 102B 

2000/5.5/0.0 

16.50 

G1 229A 

3700/4.0/-0.5 

7.44 

G1 229A 

3700/3.5/-0.7 

7.44 

G1 229A 

3600/3.5/-1.0 

7.47 

G1 229A 

3500/4.5/0.0 

9.87 

G1 229A 

3400/5.5/0.0 

10.08 

G1 229B 

1000/3.5/-0.5 



Table 5. Evolutionary Parameters 


Name 

L/L0 

[m/H] 

Teft K 

log(g) 

Mass M0 

Age Gyr 

LHS 102A 

-2.32 

0.0 

3200 

5.0 

0.13—0.19 

0.08—10.0 

LHS 102B 

-4.05 

0.0 

1900 

6.0 

0.065—0.074 

1—10 

G1 229A 

-1.29 

-0.5 

3700 

4.0 

0.30—0.45 

0.016—0.045 

G1 229B 

-5.21 

-0.5 

1000 

3.5 

>0.007 
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2.0x10 

1.5x10 

1.0x10 

5.0x10 


2.0x10 

1.5x10 

1.0x10 

5.0x10 


LHS 102A model:32-5.0-0.0-2.0 



Wavelength [micron] 


model;32-4.5-0.0-2.0 



Wavelength [micron] 


2.0x10 


1.5x10 


1.0x10 


5.0x10 


model:33-4.5-0.0-2.0 



Wavelength [mlaron] 


Figure 1. Best fit to LHS102A using our current model grids. The parameters of the models (dotted lines) are Tgff = 3200 K, 
log(^) = 5.0, a. = 2.0 (top panel); = 3200 K, log(p') = 4.5, o = 2.0 (bottom left panel); Teff = 3300 K, log(^) = 4.5, cx = 2.0 
(bottom right panel). All models shown have solar abundances. 
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Wavelength [Angstroem] 



Figure 3. Fit to higher resolution spectrum of LHS 102A using best fit high-resolution spectrum plus a number of close-by 
models. The parameters of the models (dotted lines) are Teff = 3400 K, log(^) = 5.0, [M/H] = 0.0 (top panel); Teff = 3300 K, 
log(^) = 5.0, [M/H] = 0.0 (middle panel); Teff = 3500 K, log(p') = 5.0, [M/H] = —0.5 (bottom panel). All models shown have 
mixing lengths a. = 1.5. 
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2.5x10 

2.0x10 

1.5x1 0 

1 .0x10 

5.0x1 0 


LHS 102B model:1 7-6.0-0.0-cond 



Wavelength [micron] 


2.5x10 

2.0x10 

1.5x10 

W" 1.0x10 

5.0x10 

-5.0x1 0 


model: 16-6.0-0.0-cond 



Wavelength [micron] 


2.5x10 

2.0x10 

1.5x10 

W" 1.0x10 

5.0x10 

-5.0x10 


model;18-6.0-0.5-cond 



Wavelength [mlaron] 


Figure 5. Best fits to LHS102B using “AMES-cond” models. The parameters of the models (dotted lines) are Teff = 1700 K, 
log(5) = 6.0, [M/H] = 0.0 (top panel); Tgff = 1600 K, log(g') = 6.0, [M/H] = 0.0 (bottom left panel); Tgff = 1800 K, 
log(^) = 5.0, [M/H] = —0.5 (bottom right panel). 
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5.0x10 

4.0x1 0 

5.0x10 

2.0x10 

1.0x10 


GL 229A model:37-4.0-0.5-2.0 



Wavelength [micron] 


5.0x10 

4.0x1 0 

3.0x1 0 

2.0x10 

1.0x10 



Wavelength [micron] 


5.0x10 

4.0x1 0 

3.0x10 

2.0x10 

1.0x10 
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Figure 7. Best fit to G1229A using the newest model grids. The parameters of the models (dotted lines) are Teff = 3700 K, 
log(5) = 4.0, [M/H] = —0.5, a. = 2.0 (top panel); Tgff = 3700 K, log(5') = 4.0, [M/H] = —0.5, a. = 1.5 (bottom left panel); 
Teff = 3600 K, log(5r) = 4.5, [M/H] = —0.5, a = 1.0 (bottom right panel). 








xIO 


Gl 229A (Hi-ResJ; 57-4.5-1.0-2.0 
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Wavelength [Angstroem] 



Figure 9. Fit to higher resolution spectrum of G1 229A using best fit high-res spectra plus a number of close-by models. The 
parameters of the models (dotted lines) are Tgff = 3700K, log(5 ') = 4.5, [M/H] = —1.0, ck = 2.0 (top panel); Tgff = 3600 K, 
log(5) = 4.0, [M/H] = —0.5, a. = 1.5 (middle panel); Tgff = 3700 K, log(5') = 4.5, [M/H] = —0.0, a. = 2.0 (bottom panel). 
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1.0x10 


5.0x10 



Wavelength [Angstroem] 




Figure 11. Fit to higher resolution spectrum of G1 229B using best fit high-resolution spectrum plus a number of close-by 
models. The parameters of the models (dotted lines) are Teff = 900 K, log(^) = 3.5, [M/H] = —0.5 (top panel); Tgff = 900 K, 
log(^) = 4.0, [M/H] = —0.5 (middle panel); Teff = 1000 K, log( 5 r) = 3.5, [M/H] = —0.0 (bottom panel). 
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1.0x10 
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Figure 13. Fit to higher resolution spectrum of G1 229B using best fit high-resolution spectrum plus a number of close-by 
models. The parameters of the models (dotted lines) are Teff = 900 K, log(^) = 3.5, [M/H] = —0.5 (top panel); Tgff = 900 K, 
log(^) = 4.0, [M/H] = —0.5 (middle panel); Teff = 1000 K, log( 5 r) = 3.5, [M/H] = —0.0 (bottom panel). 


































































































